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We propose a full-vectorial numerical study of far-field radiation patterns in coherent 
anti-Stokes Raman scattering (CARS) microscopy. We emphasis the particular role of the 
Raman depolarisation ratio of the observed medium and show how it modifies the radiation 
pattern of thin objects. 
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1 Introduction 

Predicted in 1965 by Maker and Therune [1], coherent anti-Stokes Raman scattering 
(CARS) was first used in microscopy as a contrast mechanism in 1982 by Duncan et al. [2]. 
It has revealed to be a powerful non-invasive tool to probe the chemical composition of 
microscopic objects. Since its renewal in 1999 under collinear configuration [3], it has been 
considered as valuable to image biological samples. CARS is a third-order nonlinear effect 
in which two beams at optical pulsations uj p and u s mix in the studied medium to give 
rise to the so-called anti-Stokes radiation at pulsation uj as = 2uj p — uj s . When to p — u s 
equals a vibrational pulsation of the medium, the anti-Stokes signal is enhanced and vi- 
brational contrast is thus generated. CARS is classically described by the third order 
nonlinear tensor ■ As a coherent process, CARS efficient generation is very sensitive 
to the phase-mismatch Ak = 2k p — k s . Several geometries have been proposed to relax 

* david . gachet ©fresnel . fr 

t nicolas . sandeau@fresnel.fr 

* herve . rigneault ©fresnel . fr 



1 



the phase-matching condition, among them collinear geometry under tight focusing condi- 
tion [4] and BOXCARS geometry [5], the former being the most implemented nowadays. 
To provide good axial and lateral resolutions, microscope objectives with high numerical 
apertures (NA) are commonly used [6,7,8]. 

As a coherent process, CARS generation is very sensitive to both size and shape of 
imaged objects. Image formation process in CARS microscopy is thus narrowly bound 
to the far-field CARS radiation pattern of studied samples, so that images are not the 
simple convolution of the object with the microscope point spread function as in fluores- 
cence microscopy. Moreover, under tight focusing condition, the commonly used paraxial 
approximation breaks. Hashimoto and al. [9] first derived the coherent transfer function 
and the optical transfer function of a CARS microscope under scalar assumption. Based 
on the framework developed by Richards and Wolf to treat the problem of tightly focused 
beams [10], Volkmer and al. solved the problem introducing Hertzian dipoles [11] and later, 
Cheng and al. used a Green's function formalism [12]. These two studies took into account 
(i) the vectorial nature of the pump and Stokes exciting fields and (ii) both the size and 
the shape of the imaged object. However, assumptions on the exciting fields polarisation 
near the objective focus were made. In particular, their longitudinal components (along 
the optical axis) were neglected, neglecting in the same time the longitudinal component 
of the nonlinear polarisation responsible for CARS radiation. 

As we will show in this paper, the relative amplitude of this component is a function 
of (i) the way the incident beams are focused into the sample and (ii) the Raman depo- 
larisation ratio (pn) of the probed medium. In other circonstances, this ratio is found 
to play an important role in elimination of the non-specific signal in polarisation CARS 
spectroscopy [13,14] and microscopy [15]. The nonlinear polarisation acting as a source 
for the anti-Stokes CARS generation, introducing its longitudinal component potentially 
affects the far-field radiation pattern of the studied sample. 

This paper starts with some CARS basics and expressions of the induced nonlinear po- 
larisation as a function of the Raman depolarisation ratio are derived. In a second part, the 
computation method used in this paper is briefly described, acompanied by a description 
of the simulated physical situation. Then, exciting fields and nonlinear polarisations are 
computed for different focusing conditions and different values for the Raman depolarisa- 
tion ratio. Finally, the influence of these parameters on far-field CARS radiation patterns, 
for different classes of objects, are presented. 



2 Raman and CARS background 

CARS is governed by the third order nonlinear tensor x^- It is the superposition of a 
vibrational resonant term, referred as x^r an d an electronic nonresonant term denoted 
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Xnr [16]- At exact vibrational resonance, the former is a purely imaginary number while 
the latter can be considered as real [17]. In the case of an isotropic medium, x depends 
on three independent components Xxxyy, Xxyxy and Xxyyx- Typical CARS experiments 
involving only two input beams, the pump field is frequency-degenerated. As a consequence, 
the number of its independent components reduces to two so that [18] 

X% = xi% y {SiM + S ik Sjl) + X^yJuSjk (1) 

where subscripts i, j, k and / refer to cartesian coordinates x, y or z, and 5 refers to the 
Kronecker delta function. 

The link between the two components Xxxyy and Xxyyx is not straightforward and it is 
useful to connect their values with experimental data obtained with spontaneous Raman 
spectroscopy. It is well known that in spontaneous Raman spectroscopy, a depolarisation 
ratio pr can be defined. It refers to the faculty of the probed Raman line to depolarise the 
polarised excitation beam. It is defined by 



where I s (\\) and / S (-L) refer to Stokes intensity respectively polarised parallel and per- 
pendicular to the excitation polarisation, p^ is specific to both the probed Raman line and 
the excitation conditions so that it can be expressed as a function of intrinsic parameters 
of the Raman line [19] 

PR ~ 45a 2 + 4 7 2 (3) 

where a, 7 S and 7 a respectively refer to the isotropy, and the symmetric and anti- 
symmetric parts of the anisotropy of the usual Raman tensor. By analogy, a, 7 S and % 
coefficents can be defined for CARS scattering [19]. 

In the same way, a CARS depolarisation ratio Pcars is defined by [19] 

(3) (3) 
Xxyyx _ Xxyyx 

Pcars - ~ 9 (3) , (3) • W 

^Xxxyy ~r Xxyyx 

Pcars is simply related to the a, % and % coefficients by the relation 

-57a 2 + 3 7s 2 



I |2 

I Pcars I 



45a 2 + 4 7s 



(5) 
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In the case when no direct electronic absorption occurs, % equals zero and both a and 
% are real. Moreover a and j s can be safely identified to their spontaneous counterparts 
(ie a = a and j s = j s ) [20]. Therefore Pcars can be assumed to equal pr. When a equals 
zero, the Raman line is told to be depolarised and Pcars equals 0.75. In the opposite case 
(7s = 0), the Raman line is totally polarised so that Pcars equals 0. Finally, under the 
assumption of no electronic absorption from the medium, Pcars lies between and 0.75. 
In the particular case of nonresonant CARS, Xxxyy equals Xxyyx in virtue of Kleinman's 
symmetry [21], and Pcars equals 1/3. 

(3) (3) 

Expressing Xim as a function of Xxxyy and pn, it is straightforward to write 



Xijkl = Xilyy(SijSkl + SikSjl + _ SilS jk ). (6) 

PR 

The local third-order nonlinear polarisation induced, at the point r, by the pump and 
the Stokes fields E p and E s is expressed by 



P (3) (r, -u as ) = x {3 '{-u as ]uj p ,u p , -u s )E p (r,u p ) : E p (r,u p ) : E*(r,-w s ) (7) 

where u p , oj s and u as are the respective angular frequencies of the pump, Stokes and 
anti-Stokes fields, the symbol * is used for the complex conjugation and the symbol : indi- 
cates tensorial product. This nonlinear polarisation is the source of the anti-Stokes field. 

Taking into account the pump field frequency- degeneracy and omitting frequency ar- 
guments uj p , uj s and uj asj the i-th component (i=x,y,z) P/ 3 "* of the nonlinear polarisation 
p( 3 ) can be expressed as 

if } (r) = 3 X%E Pj (r)E Pk (r)^ (r) (8) 
where the subscripts j, k, I run over x, y, z. 



3 Computing method and simulated physical situa- 
tion 

We have investigated the effects of tightly focused excitation beams on CARS generation 
with a fully vectorial model. The full description of this model can be found in refer- 
ence [22]. For convenience, we briefly sum up its main features. It bases on the framework 
developped by Richards and Wolf [10] to treat cases when the paraxial approximation 
breaks. Exciting pump and Stokes beams are assumed to be gaussian and are described as 
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a superposition of plane waves that are focused through a high numerical aperture (NA) 
microscope objective. The finite size of the back aperture of the objective was also taken 
into account via a parameter (3. First proposed by Hess et al. [23], it equals the ratio of 
the back aperture radius r to the half width at half maximum (HWHM) a of the gaussian 
incident beams. The resultant electric fields E p and E s , considered as vectorial, are then 
computed in the vicinity of the focal plane. They induce dipoles in the active medium 
(ie the medium emitting CARS radiation), which orientation, phase and strenght are de- 
termined by the mean of Eq.(7). These dipoles act as sources for CARS radiation, which 
far-field radiation pattern is finally computed. 



Figure 1: Schematic of the simulated configuration, f: focal distance of the objective; -Fo: 
objective's focus; r^: radius of the back aperture of the objective; a: incident gaussian beam 
half-width at half-maximum; n: immersion liquid refractive index ; 6: angle made by extreme 
rays with the optical axis. The incident beam is linearly polarised (blue arrow) along the x-axis. 
The parameter is defined by (3 = r^/a. 

CARS generation is studied when the whole vectorial components of the electric fields 
E p and E s are taken into account in the active medium, which is assumed to be isotropic so 
that Eq.(l) holds. As shown on Figure 1, the incident pump and Stokes beam (associated 
respectively to electric fields E p and E s ) are focused in the active medium through a 
microscope objective (NA=1.2 in water). They are supposed to follow a gaussian spatial 
distribution, to propagate along the z-axis and to be linearly polarised along the x-axis 
(blue arrows on Figure 1) so that they are polarised along the x- and z-axes (the component 
along the y-axis vanishes following [10]) in the vicinity of the focal plane. The higher the 
angle 9 (see Figure 1), and hence the numerical aperture, the stronger the fields components 
along the z-axis. The induced third order nonlinear polarisations along x- and z-axes equal 



Pi 3) (r) = 3 [xSLfl. M + A.K. W E t ( r ) + 2 x2L**- i^K ( r )] (9) 
^ (3) (r) = 3 [ X f zzz El (r)El (r) + £ 2 ^ (r)< (r) + 2^ xx E Pz {v)E Px (r) E* Sx (r)] (10) 



The pump and Stokes beams wavelengths are assumed to be respectively 750 nm and 
830 nm. These values match most of the recent CARS microscopy experiments, where 




Objective 
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near infrared (NIR) lasers are used [24]. Such wavelengths prevent most active media from 
one and two-photon electronic absorption, so that identification of Pcars to pr holds and 
Eq.(6) is valid. Therefore, Pj 3 ^(r) and Pz 3 \r) can be recast under the more convenient 
form 



Given Xxxyy, P x 6> (r) andPj dj (r) are now functions of r, and Pr only. Their dependence 
on r relies on the E Px , E Pz , E Sx and E Sz field maps while pr only depends on the active 
medium. 

Eventually, throughout this paper, we assume no refractive index mismatch between 
the active medium and its environment (although it has been recently shown that refrac- 
tive index mismatch can distort CARS radiation pattern [25]) while the active medium 
dispersion is assumed to be negligible (ie n(u p ) = n(oj s ) = n{uj as ) = 1.33). 

4 Mapping the components of the nonlinear polarisa- 
tion 

As shown by Eqs.(ll) and (12), the polarisations Px 3 \r) and Pj 3 ^(r) tightly depend on (i) 
the spatial distribution of the fields components E Px , E Pz , E Sx and E Sz and (ii) the depo- 
larisation ratio pr. We will be, first, interested in the behaviour of the exciting beams near 
the focal plane. Then, we will describe the induced nonlinear polarisation as a function of 
the focusing conditions and the depolarisation ratio pr of the active medium. 



4.1 Exciting fields 

Given the high numerical aperture of the objective, the exciting beams are diffraction- 
limited in the vicinity of the focal plane, following an Airy pattern. The spatial distri- 
bution of the fields components E Px , E Pz , E Sx and E Sz only varies with the parameter (3. 
As schematized on Figure 1, for any value of j3, the depolarisation of the incident electric 
fields is maximal in the (xz)-plane and null in the (yz)-plane, so that both the pump and 
Stokes exciting fields along the z-component are stronger in the former plane than in the 
latter. Consequently, for clarity, the study of the exciting fields (and a fortiori the induced 
nonlinear polarisation) will be restricted to the (x2;)-plane (although it is computed every- 



P^(r,p R ) 




where) . 
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Following Ref. [10], the exciting fields are, excepted in some particular planes, ellipti- 
cally polarised near the focus. However, throughout this part, we will only be interested 
in the amplitudes of the E Px , E Pz , E Sx and E Sz components of the fields, as computed in 
Ref. [22]. The spatial distribution of the E Px and E Pz components amplitude of the pump 
exciting field, near the focal plane, for different values of (3 are depicted in Figure 2. The E Pz 
component is rigorously null in the (yz)-plane, in agreement with Richards and Wolf [10]. 
Moreover, it is antisymmetric with respect to the focal point (E Pz (r) = E*(—r)). When 
(3 varies from 0.1 (a,d) to 1 (c,f), the tightness of the focusing decreases so that both the 
lateral and axial dimension of the focal volume increases (this effect is prevailing along the 
axial dimension). In parallel, the axial component E Pz gets lower. Such a behaviour can 
be easily explained by the "filling" of the objective back aperture by the incident pump 
beam. When (3 equals 0.1, the incident pump beam overfills the microscope back aperture, 
so that it can be considered as a plane wave. On the contrary, when (3 equals 1, it underfills 
the objective back aperture. f3 equalling 0.5 (b,e) can be considered as realistic when the 
incident beam matches the objective back aperture, what is fulfilled in most experiments. 
The same conclusions can be drawn to the Stokes exciting field. 




x position (urn) x position (fim) x position 




x position {\im) x position (urn) x position {\xm) 



Figure 2: normalised x- (a-c) and (d-f) z-component amplitude maps of the pump exciting field 
in the vicinity of the focal plane. The paprameter (3 equals 0.1 (a,d), 0.5 (b,e), and 1 (c,f). For 
each value of f3, the amplitude of each component is normalised relative to the maximum of the 
total amplitude of the pump exciting field. 
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4.2 Nonlinear polarisations 



As shown previously, the parameter (3 influences the exciting fields spatial distribution, 
and, therefore, the induced nonlinear polarisation. Furthermore, the observed medium it- 
self influences the nonlinear polarisation via the Raman depolarisation ratio pr. 

To study the influence of the observed active medium, we have considered three values 
for pr (0, 0.33 and 0.75). When pr equals 0, the Raman line is told to be totally polarised. 
Indeed, from Eqs.ll and 12, and the symbol • denoting the complex scalar product, the 
induced nonlinear polarisation is expressed by 



so that the induced nonlinear polarisation is oriented in the exciting pump field direc- 
tion. p R equalling 0.33 corresponds to a midly polarised Raman line, where Xx} yy = X^yyx- 
Eventually, in the case where pr equals 0.75, the Raman line is depolarised. Following 
Eqs.(ll) and (12), a modification of the pR value affects the respective contributions of the 
E Px , E Pz , E Sx and E Sz fields to the x- and z-components of the induced nonlinear polarisa- 
tion. To make comparison between cases where pR varies, it is necessary to normalise the 
nonlinear polarisation distribution maps for each value of pr. 

To fully understand the normalisation procedure, it is important to note that, at the 
focus, the pump and Stokes exciting fields are only polarised along the x-axis, ie E Pz and 
E Sz are rigorously null. Thus, following Eqs.(ll) and (12), the induced nonlinear polarisa- 
tion is oriented along the x-axis too, so that at this point, it does not depend on pr. 

Therefore, the expressions of the normalised components P^norm^) and Pj 3 norm(r) are 




(13) 



given by 




(14) 



A simple derivation of Eqs.(ll) and (12) gives 




(15) 




(16) 



From Eqs.(15) and (16), it can be straightforward written 
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(3), 



,apj 3) M 



dp 



R 



-I- 



dp 



R 





EsM 







;i7) 



In the vicinity of the focal plane, the E Sx amplitude being often stronger than the E Sz 
amplitude, the left side of Eq.(17) often lies between and 1. In other words, the am- 

(3) (3) 

plitude of the P x component, at a given point r, varies more quickly than the Pi one 
with p R . On Figure 3, Px norm and Pznorm are mapped in the (:rz)-plane for (3 = 0.5 and 
increasing values of pr. While P x norm does not show any significant modification when p R 
varies from (Figure 3a) to 0.75 (Figure 3c), the maximum of Pznorm decreases from 0.22 
(Figure 3d) to 0.07 (Figure 3f). 



This decay is accompanied by a deformation of the P z norm spatial distribution. Starting 
with two regular side lobes (Figure 3d), it exhibits four lobes (Figure 3e) and is finally 
butterfly-like-shaped (Figure 3f). 



(C) P H =0.75 






x position {\xm) x position {\xm) x position {\xm) 



Figure 3: normalised x- (a-c) and (d-f) ^-component amplitude maps of the nonlinear polarisation 
in the vicinity of the focal plane for (5 = 0.5. The depolarisation ratio p R equals (a,d), 0.33 
(b,e), and 0.75 (c,f). For each value of p R , the amplitude of each component is normalised relative 
to the amplitude of the nonlinear polarisation at the focus. 



To illustrate the influence of (3 on the induced nonlinear polarisation, on Figure 4 are 

(3) 

dislayed P 2 normalised amplitude profiles along the x-axis, for different values of (3 and 
Pr (the same profiles relative to the P x component are not plotted due to the weak in- 
fluence of (3 and pr). This axis corresponds to the dashed lines on Figure 3(d-f). It can 
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CL -1.0 -0.5 0.0 0.5 1.0 CL -1.0 -0.5 0.0 0.5 1.0 CL -1.0 -0.5 0.0 0.5 1.0 



x-axis position (|jm) x-axis position {[im) x-axis position (|iim) 

Figure 4: normalised amplitude profiles of the nonlinear polarisation z-component along the 
x-axis when (5 equals 0.1 (a), 0.5 (b) and 1 (c), for pn equaling (dark traces), 0.33 (blue traces) 
and 0.75 (red traces). These profiles are taken along the dashed lines in Figure 3. 

be noted that in the case where incident plane waves are focused on the sample (smallest 
value of j3), p^ has only little influence on the induced nonlinear polarisation (Figure 4a). 
However, this influence gets stronger as f3 grows (Figures 4b and 4c), although exciting 
fields intensity maps (see Figure 2) exhibits only few modification. 



5 Consequences on CARS far-field radiation patterns 

The previous study has shown the dependency of the spatial evolution of the nonlinear 
induced polarisation to both the parameter f3 and the Raman depolarisation ratio. Of 
course, as a coherent process, CARS generation is very sensitive to changes to the nonlin- 
ear polarisation. We show, in this part, how far- field radiation paterns are modified with 
Pr. As depicted in Figure 4, the nonlinear polarisation only slightly changes with pp> when 
the incident exciting beams behave as plane waves (j3 = 0.1). Therefore, the influence of 
the Raman depolarisation ratio to CARS radiation patterns has only been achieved for 
(3 equalling 0.5 (the case where f3 equals 1 seems to us quite far from usual experimental 
conditions) . 

Two kinds of objects have been investigated: thick and thin objects. The first class is 
predicted to only radiate in the forward (same direction of propagation as exicitng beams) 
direction while the second class also radiates in the epi direction (opposite direction of 
propagation to exciting beams) [12]. Of course, fine features of these patterns depend on 
the exact shape of the objects. Following Figure 3, the nonlinear polarisation is appreciable 
in the focal plane only in a lpm x lpm wide square. For this reason, we find it relevant to 
assign this transverse dimension to both objects. They are taken as parallepipeds which 
axial extent (along the z-axis) can be varied. The z-component of the nonlinear polarisa- 
tion being rigorously null at the exact focus, much smaller objects could not experience any 
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modification of their radiation pattern with changing On the contrary, larger objects 
have their emitting area limited by the size of the excitation volume. 

To simplify the problem, we will first neglect the nonresonant part of the tensor. 
In a second part, we will take into account this contribution and show how it modifies the 
radiation patterns. 



5.1 Purely resonant sample 

Neglecting the nonresonant part of X i we start with a 500nm thick object. Its forward 
radiation patterns, in the reciprocal space (k x ,k y ), are displayed on Figure 5, for p R lying 
between (a) and 0.75 (c). In this case, the Raman depolarisation ratio has little effect on 
the far-field CARS radiation pattern. The only feature to be noted is the slightly decreas- 
ing divergence of the anti-Stokes beam in the (yzj-plane as pn increases. 




Figure 5: Far-field radiation pattern in the reciprocal space (k x ,k y ) from a thick-shaped object 
(lpm x lpm, 500ra?n thick) centered in the (xy)-plane when pn equals (a,d), 0.33 (b,e) and 
0.75 (c,f), for f3 = 0.5. Each diagram corresponds to forward-emitted signal, ko = 2tt/X. Each 
radiation pattern is normalised. 

The forward and epi radiation patterns of a slice-shaped object, for various values of 
Pr, in the (xz)-plane and in the reciprocal space (k x ,k y ) are shown on Figures 6 and 7. 
The object is now an infinitely thin slice of dipoles. It is morphologically identical to bi- 
ological membranes found in cells. Following Figure 6, the radiation pattern tends to be 
symmmetrical when pn approaches 0.75 (Figure 6c). Reminding previous results obtained 
for the induced nonlinear polarisation (Figure 3), the observed symmetry conveys the x- 
orientation of the dipoles. Further information is drawn from Figure 7. First, as in the 
case of the thick medium, a very slight change in the forwardly-emitted anti-Stokes beam 
divergence is observed (Figure 7a-c). Then, a more important change in the divergence of 
the epi-emitted anti-Stokes beam occurs in the (xy)-p\a,ne (Figure 7d-f). Contrarly to the 
case of the forwardly-emitted beam, the divergence increases with p R . Figure 8 displays the 
ratio of forward to epi-collected intensity as a function of the forward-collection NA (noted 
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F/E ratio). The epi-collection NA is supposed to be constant and to equal the excitation 
NA, ie 1.2, and the forward-collected intensity is normalised relative to the epi-collected 
intensity. Naturally, for any value of pr, the F/E ratio is an increasing function of the 
forward-collection NA. Moreover, for any value of the forward-collection NA, the higher 
the pr value, the smaller the F/E ratio. A further analysis shows a relative variation of 
the F/E ratio lying between 26% and 36% when pr varies from to 0.75. It lies around 
26% for low NA (typically less than 0.3). For commonly used 0.5 NA condensors, it equals 
28% and when the collection is insured by another high NA objective (1.2 in water for 
example), it reaches 35%. 




x position (|jm) x position (|iim) x position (|im) 

Figure 6: Far-field radiation pattern in the (x2;)-plane from a slice-shaped object (lpm x lpm) 
located in the (xy)-plane when pr equals (a), 0.33 (b) and 0.75 (c), for f3 = 0.5. Each radiation 
pattern is normalised. 



5.2 Influence of the nonresonant part 

The previous investigation presented the advantage to physically highlight the main mod- 
ifications of the raditaion pattern with the parameter pr. However, it failed to depict a 
realistic CARS experiment as it neglected the nonresonant part Xnr °f the nonlinear ten- 
sor. To take it into account, we have considered the previous "thin" object, located in the 
(:n/)-plane, whose third order nonlinear polarisation is now given by 

P (3) (r) = a ■ P (3) (r, p R = 1/3) + (1 - a) ■ P (3) (r, p R = 0) ■ exp(* tt/2) (18) 

where P^ 3 ^ (r, pr) is defined by its x- and z-components in Eqs.ll and 12 and a is a 
weighting coefficient. 

p( 3 )(r, pn = 1/3) stands for the nonresonant contribution while P^ 3 ^(r, pr = 0) stands 
for the resonant contribution (note the ir/2 dephasing, with respect to the nonresonant 
contribution, at resonance). Pr = was chosen for the resonant contribution since it ex- 
hibits the strongest F/E asymmetry. 



12 




Figure 7: Far-field radiation pattern in the reciprocal space (k x ,k y ) from a slice-shaped object 
(lpm x lpm) located in the (xy)-plane when pr equals (a,d), 0.33 (b,e) and 0.75 (c,f), for 
(3 = 0.5. (a-c) forward-emitted signal (F); (d-f) epi-emitted signal (E). ko = 2tt/X. For each value 
of pr, each radiation pattern is normalised relative to the forward intensity at (k x = 0,k y = 0). 

Starting with only the resonant part (a = 0), the nonresonant part was increased from 
a tenth of the resonant part (a = 1/11) to twice (a = 2/3). The intensity ratii F/E, for 
these values of the nonresonant part (NR), as a function of the numerical aperture of the 
forward collection are plotted on Figure 9. As expected, the epi and forward radiation 
patterns of this object come closer to those of a "purely nonresonant object" (ie which Ra- 
man depolarisation ratio equals 1/3) with increasing contribution of the nonresonant part. 
The nonresonant part thus attenuates the slight differences observed in radiation patterns 
for changing values of the Raman depolarisation ratio. This is also true for objects with 
various shapes and Raman depolarisation ratii. 

In biological samples, where the imaged samples are surrounded by solvent such as 
water, we predict (for thin objects) that the forward radiation pattern is governed by the 
Raman depolarisation 1/3 of the nonresonant surrounding medium while the epi radiation 
pattern is driven by both the Raman depolarisation ratio of the object and the relative 
strength of its nonresonant part. 
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6 Conclusion 



Through this paper, a further investigation of far-field CARS radiation patterns under 
tight focusing conditions has been lead through a full-vectorial study. It has revealed the 
conjoined role of focusing conditions (through the parameter (3) and the Raman depolar- 
isation ratio pr of the studied medium, in addition of those, already known, of the size 
and shape of imaged objects. While the far-field radiation pattern of thick object is not 
affected by changes in value of (3 and pr, those of thin objects are slightly modified, con- 
cerning both the anti-Stokes beam divergence and the ratio of epi to forward-generated 
power. Such effects cannot be seen when neglecting the longitudinal components of the 
exciting fields. However, they might be only observable for strong Raman lines (such as 
the relative nonresonant part is weak) which is not always the case, specially when working 
with biological samples. In most cases, this study validates the treatment of the problem 
previously proposed by Cheng and al. [12]. However, in the case of thin objects, it brings 
some corrections. Such objects are encountered when imaging biological samples, cellular 
membranes being a few nanometers thick. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Numerical aperture of the forward collection 



Figure 8: Forward on epi detected intensity ratio (F/E ratio) as a function of the numerical 
aperture of the forward collection for pR equalling (black), 0.33 (red) and 0.75 (blue). For each 
value of pr, the ratio is normalised with respect to the epi-emitted signal intensity collected with 
a 1.2 NA objective. 

We have restricted the analysis to the case of isotropic media, excited with collinearly 
polarised exciting beams. Furthermore, the non-resonant surrounding solvent or matrix 
has been neglected in the computations but its effect can be easily predicted from the 
nonresonant contribution analysis (see section 5.2). The situation is far more complex 
when taking into account the anisotropy of studied media as well as the possible ellipticity 
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Figure 9: F/E intensity ratio as a function of the numerical aperture of the forward collection 
for pa = 0. The relative nonresonant part of the nonlinear tensor varies from (black) to 2 
(orange). In red, the case of a thin object with pn which value equals 1/3. For each value of 
the nonresonant part NR, the ratio is normalised with respect to the epi-emitted signal intensity 
collected with a 1.2 NA objective. 

of the exciting beams polarisations. It can be, of course, modelled, following the same 
electromagnetic treatment. For the case of electronically resonant CARS [26], the Raman 
depolarisation ratio no longer lies between and 0.75 and varies on a larger range [20], so 
that the situation must be reexamined very carefully. 
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